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ABSTRACT 


The precision with which the focal depth may be 
determined using Love/Rayleigh wave spectral ratios 
depends on the accuracy of the models for Earth struc- 
ture and for source mechanism used in the focal depth 
calculations, Estimates of the precision of the focal 
depth determination are obtained using the partial 
derivatives of Love/Rayleigh spectral ratios with 
respect to the parameters: focal depth, shear velocity, 
dip angle, and slip angle. We find that errors due to 
imprecise knowledge of any of these parameters can be 
important in practice. 
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INTRODUCTION 


Accurate determination of focal depth trom body 
wave information is usually not possible for shallow- 
focus earthquakes because of the difficulty in identi- 
fying depth phases such as pP. Recently the problem of 
determining the focal depth of an earthquake has been 
approached using surface wave information such as 
Spectra and spectral ratios (Keilis-Borok and Yanovskaya, 
1962; Tsai, 1969; Harkrider, 1970; Tsai and Aki, 1970a, 
1970b, 1970c; and Canitez and Tokséz, 1971). Surface 
wave spectra are tunctions not only of the focal depth, 
but also of the earthquake source mechanism and the 
physical parameters of the earth, Therefore th: source 
parameters and the Earth model must be derived from 
other information before the focal depth can be 
estimated using surface wave Spectra. The purpose of 
the present study is to determine how sensitive the 
ratio of Love to Rayleigh wave spectral amplitudes is 
to changes in the source depth, in the shear velocity 
within the Earth near the source, and in the source 
orientation, The accuracy of the focal depth determined 
from surface wave spectral ratios may then be specified 
in terms of the accuracy of the source mechanism and of 
the Earth model, 


The theory is now available for the Calculation of 
the far-field displacement associated with the propa- 
gation of surface waves generated by mathematical rep- 
resentations of earthquake source models within realistic 
Earth structures, The theory for calculating surface 
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wave dispersion in a multilayered medium was developed 
initially by Haskell (1953), and was expanded by the 
work of Harkrider and Anderson (1962), Rosenbaum (1964), 
Thrower (1965), Dunkin (1965) and Saito (1967). Mathe- 
matical representations of source models for earth- 
quakes were derived through the work of Yanovskaya 
(1958), Knopoff and Gilbert (1960), Ben-Menahem (1961), 
Maruyama (1963), Haskell (1963, 1964), Burridge and 
Knopoff (1964) and Haskell (1966). The relations yield- 
ing the far-field surface wave displacements due to 
source models located in specified realistic Earth 
structures were presented by Harkrider (1964) and 
Ben-Menahem and Harkrider (1964), and have been expanded 
upon in the work of Harkrider and Anderson (1966), 
Saito (1967), and Harkrider (1970). Using the theory 
developed in these studies, spectral ratios of surface 
waves are determined in the present study for a series 
of source mechanisms, focal depths, and Earth struc- 
tures. An estimate is then made of the precision of 
focal depth determinations from surface wave spectral 
ratios, 


THEORY 


Using the notation of Harkrider (1970), the far- 
field expressions for Rayleigh and Love wave spectral 
amplitudes may be written 
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where S is the spectral source function, m = 0 for a 
point force source and m = 1 for a couple of double- 
couple source, kp = w/c, for Rayleigh waves, where Cp 
is the Rayleigh wave phase velocity, ky = w/C) for Love 
waves, where Cy is the Love wave phase velocity, h is 
the source depth, 6 is the azimuth from source to 
station, r is the epicentral distance, and YR and Mi. 
are the Rayleigh and Love wave attenuation coefficients 
respectively. The terms ER and ED in equation (1) are 
given by 


gE = €. Ae.’ (2) 
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where Ey is the Rayleigh wave ellipticity 
eX 
ey = - 0,78) (3) 


and Ap and Ay are the Rayleigh and Love amplitude 
response due toa vertical point force at the surface, 
The radiation pattern function x(@,h) given in equation 
(1) is 


x(@,h) = d, sae | (d) Sin 6 + d, cos 8) + d, Sin 26 


(4) 
+ dy cos 26 
where the coefficients d; are defined in Table I, 
Using equations (1) and (2), the Love/Rayleigh 
amplitude ratio may be written: 
' m “1/2 2Yi, Fr 
m x -172 “Yor 
R kp R EG Ap kp e ‘R 


assuming the spectral source function S to be approxi- 
mately the same for Rayleigh and Love waves, 


To determine the sensitivity of the ratio Love/ 
Rayleigh to focal depth h as compared to its sensitivity 


to changes in the shear velocity 8 of the Earth and in 
the dip angle 6 and slip angle A of the fault, the 
following ratios of partial derivatives may be computed: 


DR(h, 8B) 


DR(h, 6) 


(6) 


DR(h,A) 


From equations (5) and (6), we have: 
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ar | fr 
DR(h,é) = » and DR(h,A) = Tan ° (7) 
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These ratios of partial derivatives are independent of 
the attenuation coefficients Yp and Y; and are depen- 
dent on frequency w and azimuth 6, 


CALCULATIONS FOR THE BASIN AND RANGE 
VELOCITY STRUCTURE 


The Basin and Range province was chosen as a source 
model for study of Love/Rayleigh spectral ratios both 
because the earthquakes occurring within this province 
are typically shallow-focus and because of the avail- 
ability of a large amount of seismic data for this 
region making possible tlhe construction of a velocity 
model which may be more accurate than that obtainable 
for most seismically active areas. The P velocity 
structure chosen to represent the Basin and Range 
source region is based on an upper mantle model (CIT 111P) 
determined by Archambeau et al., (1969). The crustal 
model is similar to that determined by Warren (1969) for 
the southeastern part of the Basin and Range province, 
The P and S velocity structures as well as the density 
distribution for the Basin and Range source model are 
shown in Figure 1 for the crust and in Figure 2 for the 
upper mantle, and the parameters for the model are 
listed in Table II, The S velocity structure was computed 
from the P velocity structure by assuming a Poisson's 
ratio of 0.25, 


Two source models were considered in the present 
study: left-lateral fault with dip angle of 90° and 
slip angle of 0° and a left lateral fault with dip 
angle of 90° and slip angle of 45°, Calculations were 
made for focal depths corresponding to a source 
successively positioned in the middle of each of the 
first six layers of the Basin and Range model. The 
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partial derivatives 3J/98 were computed for the first 
six layers for both source models. This was accomplished 
by successively incrementing the shear velocity 8 of the 
first, fourth and fifth layers by A8 = 0.1 km/sec, and, 
since the velocities of the second and third layers are 
nearly equal, incrementing the shear velocities of both 
these layers together by Ag = 0.1 km/sec, Spectral 
amplitudes Uy and Up were then calculated, The partial 
derivatives of spectral ratio with respect to dip angle 
6, 33/35, were computed by varying 6 by 10° and comput- 
ing the spectral ratio for each source depth for slip 
angles i} of 0° and 45°, Partial derivatives of spectral 
ratio with respect to slip angle, aJ/9A, were calcu- 
lated from the difference in the computed spectral 
ratios for slip angles 0° and 45°, Because of the large 
differential in slip angle, the calculated derivative 

is accurate enough only for rough estimates. Depth 
partials, 3J/dh, were computed by taking differences in 
spectral ratios for sources within the jth and (i + 1) th 
layers, It is possible by means of the depth derivatives 
given by Harkrider (1970) to determine analytically 

the partials 3J/dh for changes in source depth which 
occur completely within any given layer, 


The fundamental mode phase and group velocity 
dispersion for Rayleigh and Love waves for the Basin 
and Range model are shown in Figures 3 and 4, The 
fundamental mode radiation pattern functions Xp and 
Xp, for the two sources studied are presented in 
Figures 5 through 20 for periods of 20 to 40 seconds, 
These radiation patterns show that the amount of energy 
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propagating from a source decreases with increasing 
depth for periods of 20 to 40 seconds, For the source 
with a strike of 0°, a dip angle of 90° and a slip 
angle of 0°, both the Rayleigh and Love radiation 
patterns for periods of 20 and 40 seconds have four equal 
lobes for all depths. For the source with a strike of 
0°, a dip angle of 90° and a Slip angle of 45°, only 
the Love wave radiation patterns have four equal lobes 
Over the depth range examined. Love/Rayleigh spectral 
ratios are given in Figures 21, 22 and 23 for a few 
selected source geometries, Figure 21 illustrates the 
variation in J as a function of 6. The dependence of 
J on the source parameters \ and 6 is illustrated in 
Figure 22, and Figure 23 shows the change in J due to 
change in only focal depth, From Figures 21, 22 and 23, 
the orientation of the source can be seen to have as 
great an effect on the Love/Rayleigh spectral ratio as 
does the focal depth. The period at which a maximum 
occurs in the Love/Rayleigh ratios (corresponding to a 
minimum in the Rayleigh wave amplitude) can also be 
seen to vary with source parameters such as dip and 
Slip angles as well as with focal depth, 


The partial derivative ratios DR(h,8), DR(h,6) and 
DR(h,A) given in equation (6) may be interpreted as the 
changes AB, A&, and Ad respectively which would produce 
the same change in J as a change Ah in h of 1 km, There- 
fore the values of DR(h,8), DR(h,é) and DR(h,A) given 
in Tables III, IV and V respectively may be used to 
determine how accurately the shear velocity of the 
Earth and the dip and Slip angles of the source must be 
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known in order to determine the focal depth to some 
desired accuracy, Consider, for example, that it is 
desired that the focal depth be determined correct to 
within + 5 km, Then the values of AB, Ad, and Ad which 
would produce the same Change in J as a change in h of 
5 km are given by multiplying the values of DR(h,8), 
DR(h,5) and DR(h,A) in Tables III, IV and V by a 
factor of 5. The resulting AB, AS and AN are larger 
respectively than 0,2 km/sec, 10° and 10° for sources 
in certain layers but are smaller for sources in others, 
Moreover, AB, Ad and Ad are a function of azimuth and 
if a combination of Changes in shear velocities, dip 
angle, and slip angle were made, the result could 
equal to a Ah of 5 km with AB < 0.2 km/sec, AS < 10° 
and 4d < 10° for sources within several of the layers. 


Therefore the conclusion must be that the accuracy 
of the focal depth determined from Love/Rayleigh spec- 
tral ratios is very much dependent not only on the 
velocity structure and source parameters but also on 
the layer in which the source is located and the 
azimuth of the recording station from the source. In 
Some cases, the shear velocity of the model near 
the source would need to be known to better than 
+ 0.2 km/sec and the dip and Slip angles to better 
than + 10° in order to attain an accuracy of + 5 km 
in the focal depth. 


If spectral ratios are available for several 
different values of period and several Stations, the 
depth resolution could be better. However, corsideration 
of the first or fifth column of Table III for a slip 
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angle of 0° shows that the error can be of the same 
Sign at every station and frequency; and one would not, 
therefore, expect it to be minimized by any "least 
Squares" procedure, With respect to determining depth 
by finding nulls in the snectra, column two of Table III 
Shows that the errors can have the same pattern as a 
function of azimuth at every frequency, thus leading to 
a consistent shift of the null frequency, 
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DISCUSSION 


The accuracy of the focal depth determined from 
the spectral ratio of Love wave to Rayleigh wave energy 
depends upon the accuracy of the mathematical models 
used for the Earth structure and for the source 
mechanism, Recent studies by Tsai and Aki (1970a, b, c) 
and Canitez and Toksoz (1971) have concluded that it 
is possible to determine the focal depth to an accuracy 
of a few kilometers if the source mechanism and the 
Earth structure are reasonably well known. However, in 
order to define the terms "reasonably well known", the 
partial derivative DR(h,8), DR(h,é) and DR(h,A) should 
be available for models of the Earth structure and for 
the source mechanism which are close to those being 
assumed in the calculation of the focal depth. In some 
cases it probably will not be possible to estimate the 
focal depth to within a few kilometers of the correct 
value from surface wave spectral information alone. 


We see that it is not true, as often hypothesized, 
that error due to lack of knowledge of the velocity 
structure will be reduced to negligible proportions by 
using Love/ Rayleigh spectral ratios, instead of 
Love or Rayleigh spectra values alone, to determine 
source parameters, Also the difficulty in determining 
the dip and slip angles to better than 10° may in many 
cases preclude the determination of focal depth from 
surface wave spectral information to an accuracy of a 
few kilometers. 
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Figure 5. Rayleigh radiation pattern for a double couple 
representation of a left lateral vertical fault source at 
a depth of 0.5 km and with a strike of 0° and a slip of 0°, 
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Figure 6, Rayleigh radiation pattern for a double couple 


representation of a left lateral vertical fault source at 
a depth of 5.5 km and with a strike of 0° and a slip of 0°, 
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Figure 7, Rayleigh radiation pattern for a double couple 


representation of a left lateral vertical fault source at 
a depth of 14,75 km and with a strike of 0° and a slip of 0°, 
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-1gure 8, Rayleigh radiation pattern for a double couple 


representation of a left lateral vertical fault source at 
a depth of 25.75 km and with a strike of 0° and a slip of 0°, 
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Figure 9, Rayleigh radiation pattern for a double couple 


representation of a left lateral vertical fault source at 
a depth of 34.75 km and with a strike of 0° and a slip of 0°. 
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Figure 10, Rayleigh radiation 
representation of a left later 
a depth of 42,75 km and with a 


pattern for a double couple 
al vertical fault source at 
strike of 0° and a slip of 10°, 
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Figure 11, pelagaen radiation pattern for a double couple 


representation of a left lateral vertical fault source at 
a depth of 0.5 km and with a strike of 0° and a slip of 45°, 
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Figure 12, Rayleigh radiation pattern for a double couple 
representation of a left lateral vertical fault source at 
a depth of 5.5 km and with a strike of 0° and a Slip of 45°, 
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Figure 13, Rayleigh radiation pattern for a double couple 
representation of a left lateral vertical fault source at 
a depth of 14,75 km and with a strike of 0° and a Slip of 45°, 
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Figure 14, Rayleigh radiation pattern for a double couple 


representation of a left lateral vertical fault source at a 
depth of 25,75 km and with a strike of 0° and a slip of 45°, 
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Figure 15, Rayleigh radiation pattern for a double couple 
representation of a left lateral vertical fault source at a 


depth of 34,75 km and with a strike of 0° and a slip of 45°, 
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Figure 16. Rayleigh radiation pattern for a double couple 
representation of a left lateral vertical fault source at a 
depth of 42,75 km and with a strike of 0° and a slip of 45°, 
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Figure 17, Love radiation pattern for a double couple 
representation of a left lateral vertical fault source at a 
depth of 0.5 km and with a strike of 0° and a Slip of 0°, 
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Figure 18, Love radiation pattern for a double couple 
representation of a left lateral vertical fault source at 
a depth of 42.75 km and with a strike of 0° and a slip 
of 0°. 
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Figure 19, Love radiation pattern for a double couple 
representation of a left lateral vertical fault source at 
a depth of 0.5 km and with a strike of 0° and a slip 

of 45°, 
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Figure 20, Love radiation pattern for a double couple 
representation of a left lateral vertical fault source at 
a depth of 42,75 km and with a strike of 0° and a slip 
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Figure 21. Love/Rayleigh amplitude ratio as a function of 
period and azimuth, 
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Figure 22, Love/Rayleigh amplitude ratio as a function of 
period, slip angle and dip angle. 
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